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ABSTRACT 


The preparation of the analytical approach and of the related softtvare to be used 
in the inversion of the differential and rotating doppler data that the ionospheric 
e:q)eriment of the ApoUo-Sc^uz Test Project (ASTP) will obtain has been recently 
completed. These data will be collected in space-to-space paths (between the ASTP 
Docking Module (DM) and the Apollo Command Service Module that will orbit both at 
a height of 221 km, at a relative distance of 300 km) and in space-to-ground paths 
(between the DM and ground). The dealer links will c^erate at 162 and 324 MHz and 
will have an accuracy better than 3 mHz over 10-sec Integration time. The ASTP 
Mission is scheduled for summer 1975. While waiting for the actual experiment data, 
the inversion approach has been tested with the "dummy" data obtained with a com- 
puter simulation. We found that a measui*ement accuracy of 1 to 10% in the value of 
the horizontal electron density gradient at 221-km altitude can be achieved, in space- 
to-space paths. For space-to-ground paths noa»' the orbital plane, we have identified 
and compensated for possible effects of the horizontal gradients on the received 
differential d(^)pler shifts. In the two cases we have considered, it was possible to 
reduce the gradient-associated errors in the inversion that leads to the columnar 
electron content by approximately one-half. Accuracies of 5 to 10% in columnar 
electron content are achievable, with this gradient-compensation technique. 
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DOPPLER MEASUllEMENTS OF THE IONOSPHERE ON THE OCCASION 
OF THE APOLLO-SOYUZ TEST PROJECT 

* 

Part II; Inversion of Differential and Rotating Dqjpler Shifts 
Ray H. Gay and Mario D. Gross! 

1. INTRODUCTION 

The forthcoming Apollo-Soyuz Test Project (ASTP) gives us an opportunity to collect 
valuable observations cf the columnar electron content in the ionosphere at an alti- 
tude of 221 km and below. In addition, doppler measurements between the Apollo- 
S<^z docking module (DM) and tlie command service module (CSM) make it possible 
to obtain relatively continuous signatures of the electron-content gradient in the space- 
craft's orbit over a 24-hour period. After determining the aloft gradient, we can then 
consider the effect of electron-density horizontal gradients chi the DM -to-ground differ- 
ential doppler shift. Although the t<?tal gradient contribution is along the ray trajectory, 
accurate measurements of the gradient at the spacecraft altitude may be fitted to exist- 
ing ionospheric models in order to improve the descr?pticm of gradient integrals along 
the spacecraft-to-ground path. 

The detection of mass concentrations in the earth from measurements of small 
accelerations between the DM and the CSM is dependent on an accurate determination 
of the ionospheric doppler shift in the DM/CSM radio link. A 1-Hz differential doppler 
shift Is possible in the presence of large gradients and corresponds roughly to an 
apparent relative velocity of 1/3 m sec~^ between the receiver and the transmitter. 

Such a value is 100 to 1000 times as lai’ce as the expected mascon doppler signature. 

In nighttime regions, where the horizontal gradients are small, the gravitational and 
ionospheric doppler magnitudes are comparable. The use of two oscillators operating 
coherently at different frequencies enables the dynamical part of the doppler shift to 
be isolated from the ionosplieric component, but the determination of the DM/CSM 
gradient and resulting individual frequency shifts due to the ionosphere must be quite 
accurate in order to maintain the integrity of the residual mascon -induced doppler 
shift at the 2- to 3 -mHz level. 

* 

This report co\'ers the inversion of differential and rotating doppler shifts (Part II). 
Part I (SAG Special Report No. 366) co\^ers the computer simulation of the ionospheric 
experiment results. 

This research was supported by NASA under contract NAS 9-13837 to the Smithsonian 
Astrophysical Obser\*atory. 
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2. REDUCTION OF DM/CSM DOPPLER SHIFT AND RECOVERY OF 
HORIZONTAL ELECTRON-DENSITY GRADIENTS 


First, it is informative to develop an approximate e3q>ressiC'n for the differential 
dopplor shift between the two spacecraft as a function of the electron density along the 
ray path. The negative phase -path defect, 

L = P - p , 

is the difference between the straight-line distance p and the radio length P of the ray 
path. The phase-path length is related to the refractive index p by 

B 

P = f u dS , (1) 

A 

' * 

whe'.e dS is an element of length along the curv'ilinear ray path from point A to point 
B in space. 


For the ASTP experiment, the instrument frequencies (f^ = 324 MHz, ~ 162 MHz) 

are much higher than either the electron gyrofrequency f„ or the critical or plasma 

“ 5-1 

frequency f^^ in the ionosphere. For a typical electron density of 10 el cc ard n 
magnetic-Reld strength of 0. 28 Oe, = 0. 8 MHz and f^ = 2. 8 MHz, 


Neglecting collisions and the earth* s magnetic field, we can write for the refrac- 
tive index (Kelso, 1964, Chap. 5) 


p- 3 , 

2itmr 

where e and m are the electron charge and the rest mass, respectively, and N is the 
electron density in el oc~^. In mks units, we have 


jj 1 


40.3 N 


f 


where f is in Hertz, 


2 



( 2 ) 


Tlie expression for the total doppler shift between the transmitter A and a 
remo^’ed point B is 

A 

When the range p differs insignificantly from tlie arc length/ dS, the two quantities 
dp/dt and (d/dt) / dS can be interchanged without serious errors being introduced. 

In fact, tlie dual-^quency technique removes these terms from consideration in first- 
order expressions for the differential frequency shift. The differential doppler shift 
obtained from the coherent emission at two frequencies f ^ and fg is 


64> 




(3) 


and if we allow the same integration path for the two frequencies in equation (2), we 
can write 


B 




2\ 




(4) 


where the quantity I N dS is the total electron content per unit area between A and B, 
The total time derivative of J N dS, in general, is a complicated expression involving 
the motion of points A and B and various partial derivatives of N; but when the doppler 
link is between two satellites in coplanar circular orbits, the expression is 


DM CSM 

1/ NdS-N(CSM)v^.g^-N®M)Vj,„* J fdS 

CSM * DM 


(5) 


where v^gj^ and are the velocity components of the CSM and the DM in the direc- 
tion of propagation at their respective end points of the ray path. If the X axis has its 
origin at the CSM and passes through the DM and if the two spacecraft are in nearly 
coincident circular orbits, then equation {5) becomes 


^ J N dX = IN(CSM) - N(DM)] Vq + dS , 


(5a) 
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where is the common orbitnl speed. The first term in equation (5a) can be written 

AN 

IN<CSM) - N(DM)1 = Vjj ~ AX , 

80 that equation (5) is finally 

DM CSM 

1/ / f dx , 

CSM DM 

where W/8X is the average gratUent along the path coraiecting the two spacecraft, and 
SN/0t is the temporal variation of electron density in the path. If we neglect or can 
otherwise account for the 9N/3t term, tlien equation (4) gives a direct measurement of 
the electron-density gradient at the ASTP altitude averaged over the DM/CSM separa- 
tion arc. 

Thus, tlie averj^e horizontal gradient involved is given by 

® 40.3 (4-q) ''0^ • " 

In order to obtain a global view of tlwj smoothed 64> in the DM/CSM path, we com- 
puted a series of phase-path lengths for each frequency at points along the orbit sepa- 
rated by If)® in longitude (Figure 1). The ionospheric electron-density model for sum- 
mer 1975 (Odom, 1975) was used. The temporal effect of 9N/9t and the earth's magnetic 
field were excluded from the simulation. The averaged 6$ between successive points 
along the orbit corresponding to 1200 LT, 1300 LT, . . . was approximated at the mid- 
points of each longitude segment by simply taking 6$ = (A<t/A long) (A long/At). These 
1-hour intervals in geocentric longitude correspond to time intervals of 3. 45 to 4, 93 
min along the orbit (whose period is ~88 min), so they represent avsrages of many 
individual lO-sec observations of the expected 6^. 

The e3q>ected DM/CSM differential doppler shift for an entire orbit iS shown in 
Figure 2. The steepest gradients occur between 1030 and 1230 LI’ and result in a 


(5b) 
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Figure 1. Simulated ASTP orbit geometry. 
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Figure 2. Expected differential doppler shift in the DM-to-CSM path for an entire orbit. 



peak of about 300 to 400 mllz in 64>, in contrast to the 2400 to 0000 LT region, wliero 
64» is almost nogligiblo. Overall, the horizontal gradients at the 221-km ASTP alti- 
tude vary by a factor of 100 between night and day. 

In Figure 0, we show the horizontal gradient DN/DO (in el cc**^ rad~^i at the 221-km 
orbital height In the direction of Increasing (eastward) orbital longitude 0, TJie solid 
curve is plotted through the averaged values AN /AO at the midpoints of each 1-hour 
longitude intor\’al (1230 LT, 1330 LT, . . . ). The broken curve is fitted to the expected 
average values measured by the DM/CSM linl< o\'er the same 15* longitude Inten als. 
Positive values indicate increasing electron density in the direction of orbital motion. 

The e?ctreme diurnal variation in ON/flX along the orbit could reflect on the gravi- 
tational part of the experiment, in that the maximum error of tlie Ionospheric-related 
frequency shift in the individual frequencies \ s!!ies from about 2 to 200^’(;. Figure 4 
indicates the allowable error in the determination of W/DO as a function of gradient 

intensity for two values of maximum absolute error n In the received frequency of tlie 

— — 5 -1 -1 

324-MHz signal. For example, we see that for a value of 9N/90 - 6 X 10 el cc rad , 

we must have less than a 2, 2‘<1 error in the observed value of ON/aX (or 9N/0O) in orxler 

that the error in be less than 2 mHz. The nighttime logical is the most attractive 

for detecting gravity anomalies, as least as far as the ionospheric effect is concerned. 

In daylight and transition regions (evening, morning), the required accuracy approaches 

1% or less, su^esting tliat the likelihood of a small error in the ionospheric doppler 

correction masquerading as a gravity anomaly in the DM/CSM dealer record is greatly 

increased. 

The preceding observation points out the absolute necessity of employing a phase- 
coherent dual-frequency linl< to measure the ionospheric frequency shift between space- 
craft when accurate doppJer range-rate meusui*ements ai^ desired. An a priori model 
of the ionosphere simply cannot be expected to embrace the e\*ident accuracy require- 
ments of the ASTP doppler-tracking experiment. To a high degree of accuracy, the 
individual correcticais to the received frequencies f^ and fg can be calculated from 
equations (3) and (4) to be 



ACTUAL AVERAGED AN/A© 
DM/CSM LINK MEASUREMENl 
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Figure 3. The horizontal clectron-dwjaity gradient 9K/89, in el oc rad , at the 221-km orbital height. 
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Figure 4 




and 




In our case = 324 MH; . fg = 162 MHz), we have 

• 9 • 

and 


A seccmd simulation was done in order to determine in two instances just how 
well a typical averaged observation of 6# in the DM/CSM link would represent the 
ac Aial gradient at the midpoint of that portion of the orbit traversed during the actual 
integration tfn>« of the instrumentation on board the CSM. The geometry of the situation 
is seen in Figure 5, where we show the DM and ?.i o CSM separated by a constant dis- 
tance of 300 Ion as the baseline moves through 70 km in the orbital path. The 70 km 
corresponds to approximately 9 sec of real time, a figure that is close to the 10-sec 
receiver/recorder integration time. The baseline midpoint was taken to correspond 
to 0500 and 0900 LT in order to include both small and large horizontal gradients of 
electrcm densi y. By ray liacing, the averaged differential phase measurements were 
simulated and an average differential doppler shift, 


9 sec 


= =-L 

At At 


J 


64> dt 


9 


was calculated. Substituting 6^ into equation (6) then gives the average value of 9N/aX 
over the moving DM/CSM baseline. 


In Table 1, we compare values of with those of 8N/8K at the midpoints 

(0500 and 0900 LT) of the moving baseline and note the associated relative error in 
3N/8K, In both cases, the error is well below the maxim am indicated in Figure 4. 
This is encouraging, but we must also consider the possible effects of time dependence 
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508-054 



Figure 5* Geometry of Individual lO-sec dofipler-shift observations. 



in tha alaotron danstty (dN/dt ^ 0) and of large aacond-ordar darivativea d^N/£D^, 
which would dialocata tha mean pt^nt (where 1ff/8K - aN/8K) from the aaaumed mid 
point. 


Table 1. DM/C8M gradient error. 



W/ax 

0 N/aX 

1 

Errpr in afy/ax 

Limit 

Time 

(el CO 

km ) 

<%) 

(%) 

Night 0500 LT 

16.5 

14.9 

11 

136 

Day 0000 LT 

819 

823 

0.5 

2.5 
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3. INVERSION OF SPACE CRAFT -TO-GROUND DIFFERENTUL AND 

ROTATING DOPPLER SHIFTS 

S, 1 Differential Doppler Shift 

In this section, we will consider the various options available in analyzing spacecraft- 
to-ground differential doppler measurements. The effects of horizontal gradients and 
nonstationarity in the ionosphere ccanplicate the problem considerably, but it is still 
possible to extract useful Information about the ionosphere when measurements are 
made continually or at short intervals over a visible arc of the satellite orbit (Hibbcrd 
and Thomas, 1959; APpert et al , . 1965; Tyagi, 1974). The effect of the ionosphere 
on spaoecraft-to-ground doppler shifts is given precisely by the following relation: 



Assumptions inherent in the above expzession enter through equation (2). 

Considering the electron density to be a function of coordinates and time, 

N = N(x,y,z,t) , 

leads to an expression Involving those four quantities related to the differential dcpp'ter 
shift: 



following the notation of Al'perl (1965). In this equation, the satellite-bome trans- 
mitter is assumed to be moving with respect to an earth-fixed observing station through 
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a nonhomogeneous Ionosphere, where the electron density at any point alwig the ray 
path may vary with time. The cartesian coordinates (x,y, z) of a point on the ray 
trajectory are given in an earth-fixed system centered at the ground station. The 
(x, z) plane contains the station, the instantaneous satellite positicm, and the center 
of the earth. The angle 4> refers to the zenith angle of the ray trajectory, and the sub- 
script c den<ftes values at the satellite position (x , y , z ). In Figures 6 a arid 6 b, the 
relationship lietween the various quantities can be more easily visualized. It Is 
important to ncte that in the polar -coordinate system shown, the distance R 0 (or RX) 
frrnn the origin is always positive, as is the zenith angle 4 >. As a consequence, direc- 
tional derivatives of electron density 8N/90 and 3N/®C change discontinuously at the 
pole (<])q = 0 ). 

The first term within brackets in equation ( 8 ) relates the electron density in the 
vicinity of the satellite to the vertical velocity component z and is especially important 
in near-vertical rocket flights, where the vertical profile of N(z) can be easily mapped 
through observations of 64. The second term forms the ccmtributicm of the motion in the 
(R, 0) plane and defines the refraction angle 64 )^ at the observing station. The term 
- [0?J/0y] y^ gives the effect of motion normal to the (R, 0) plane and defines the cor- 
responding transverse refraction ar^le The final term, - J 3N/9t dS, accounts 
for the fact that the electron density along the ray path may be changing independently 
Mdth time. Normally, it is quite small (Al'pert, 1965) compared to the other terms 
and can be neglected unless a large-scale ionospheric disturbance occurs in the space- 
craft-to-ground path. 

The integrals and [9N/9x] are taken along the ray path F, which, to first 
order, can be considered to coincide with the geometric line of sight. They represent 
the mean electron density and the mean horizontal gradient in electron density between 
the satellite and the ground station. In terms of spherical coordinates, these quanti- 
ties are expressed as follows (R^ = earth's radius); 

J NR (R cos cj>) ^ dR 
[ 5 ^ 1 _ __________ _ ^ 

(Rq cos cj>p) - (R^ cos 4 )^) 
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( 10 ) 



1 f ^ 
Bin4>Q J 9x 


(Rq cos - (R cos ^ 

(Rq cos - (R^j cos 4.^)"'^ 


and 


® u- fM s^ne dR 

dyj ~ dy sin 9 cos <|> ^ 

* C 


( 11 ) 


where &N/3X = (1/R) \9N/80) is the horizontal gradient at a point on the ray path. If 
a planar approximation Is made (R^ — «), these expressions simplify to 



(9a) 


M . 1 f“ 8N ^az 

L0k J z^ sin 4 >q cos (i)^^ J 3x ’ 


(10a) 


and 


= 1 M 

dyj z cos 4 iq J 0y 
0 


z dz 


(11a) 


The inverslcm of the full equation (8) embraces many questions concerning just 
how to treat the obseirvational data in order to identify the different terms and obtain 
from them any definite information on the ionospheric stnicture. The local term N 

c 

is the only one that promises to give specific local measurements of electron density; 
all the others are averages along the ray path and are difficult to relate to particular 
ionospheric regions. Special techniques can be employed in some cases to reduce 
the number of unknowns in equation (8). If the observing station is located in the plane 
of the orbit, then t9N/8y] is removed from consideraticm. A nearly circular orbit 
renders insignificant, while an elliptic orbit will enable measurements to be made 
when r^ + (z^/cos (|)q) = 0, thereby removing and laR/ax]. This last condition 
has been exploited by Al’pert and Sinel'nikov (1965) to measure tBiVfiyI and N . 

C 
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3, 2 Rotating Doppler Shift 


The rotating doppler (or Faraday rotation) effect can also be used to obtain 
electron-density information if a linearly polarized antenna is utilized at the observ ing 
station. The frequency of amplitude variaticm of the circularly polarized wave as 
received is given by Al'pert (1965), where quasi-longitudinal propagation and a constant 
magnetic -field vector along the ray path are assumed. The rotating di^pler shift is 
then givon by 



Here, the minus and plus signs refer to the ordinary and extraordinary waves, 

respectively, and the signs of r , z , and y are determined in the same way as in 

c c c 

the differential dq[}pler shift. The quantities in brackets are 


ISrhI.' 


Sin Y 0 cos 


_r 

1 Rq sin <j)Q t(R^ cos c|>q)"^ - (R^ cos (R cos ^) 


NR^ dR 


(13) 


.-1 


r^l 1 r 3N cos 0 

Lax Jh " sin 4>o J ax cos cj, ^^,-1 


(14) 




sin Yo «in sin r ^ dR 


/ 


R COS ’ 


(15) 


f ^ sin 9 cos 0 
_^Jh ” ^ sin 0 ^ cos 4 ) ° 


and 


r9?i1 _ f ^ COS 0 
L9t Jjj J 9t cos 4 ) 


(16) 


(j?) 
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The following quantities are depicted in Figure 7; 

Hq = the magnetic-field vector, 

Q the angle between the magnetic-field vector and the wave normal, 

Yq <= the angle between the magnetic-field vector and the geocentric radius 
vector, 

Iq = the angle between the geometric ray of sight and the projection of onto 
the local horizontal plane. 

The upper sign in equation (15) is used when di^^/dt > 0, or whan the projectiem of IJq 
on the local horizcmtal plane is moving away from the (R,G) plane. 


As in the differential doppler relations, a planar a^Droximation clarifies the 
relationships between the quantities in equations (13) to (17) and tte physical ionospheric 
parameters. In the fi.r,t-earth case, equations (13) to (17) become 


sin Yq cos it ic 

c u 0 




cos 




- — r 

Q cos <15 q J 9x 


and 


sin Ya sin U fC 

J -j • 

n. 




COS0Q -C gjj 


as 


at 


cos 0fv /-c 

pj cos (j) 


n 


where Rq — dR dz, and the angles and 0 are constant. 


(13a) 


(14a) 


(15a) 


(16a) 


(17a) 
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toi >064 



F'lgfure 7. Geometry and notation of rotating doppler shift. 


19 



In the absence of horizontal gradients, we have for the two frequencies and a 
differential doppler shift, 6^, which can be used to solve for Jn dz and for each 
point observed during a pass. The value of N_ will be less significant because of the 

V 

near circularity of the satellite (DM) orbit. For the A8TP mission, the planned eccen- 
tricity of 0. 0002 during tlie doppler tracking phase results in an amplitude of tq^roKi- 


mately 1. 6 m sec~^ in z . Assuming N “ 10® el cc”^, this value of z would contribute 

c c c 

at most about 2 mHz to the received differential doi^ler shift. This is a negligible 

value, but an increase in a few kilometers in Hie apogee/perigee height difference could 

6 -1 

combine with an electron density of several times 10 el cc to give a nonnegligible 

value of perhaps 20 mHz owing to N z . 

c c 


.6 


-1 


If we con safely neglect the effect of N^z^, then four unknowns remain in the 
system of two equations: 

8i = f(/Ndz , Jf^d= , , 


6i„=g(jKdz , J^7.dz , 07: dz , 0^ . 


(18) 


To reduce the number of unknowns in equation (18), we laust either neglect the gradient 
integrals or make some correction to the observed quantities from a priori knowledge 
of 9N/9x and 9N/9y along the ray paths. Another condition that can be met for at least 
one observing station is that the effect of J(dN/^y) z dz can be made small in com- 
parison to the remaining terms by locating the receiver near the orbital plane. Then, 
if 9N/9t can be accounted for, the system (18) is not insoluble. 
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4. Tift: INFLUENCK OF HOIUZONTAL GRADIENTS ON THE SPACECUAFT- 
TO-GROUND DATA-INVERSION ACCLllACY 


4» 1 InverHlon In tie Prewence of Unknovsn Horl/.ontal Gradienta 


If the gradient* related terms (a?} ax) and jaS ay) are neglected !ind the time 
stationarity of the electron density N is assumed, then we can assign the total b4> to 
the effect of the path's average electron density Equation (8) then l>ecoine8 


~ { N(x,z) 
0 \ c cos<>q/ J ' » ' 


dz 


(P») 


with the planar approximation in effect. We note here that the integral (l/z^,) /n dz. 

is closely approximated by the average electron density along the straight-line path 

/* 

connecting the receiver with the transmitter, (l/p„) I N(x, z) dp, and we will use 

V irQ 

these interchangeably. Since the gradient term 1^5/ dxl can iie either negative or 
positive, depending on the geometry of the situation, we can accordingly expect the 
estimate of lS„l to ije too large or too small. However, we do not expect that gra- 
dients large enough to change the sign of 64> would exist along the spacecraft -to-ground 
path except when the satellite is almost exactly abo\e the oi)sen ing station. The 
inversion of b* to columnar electron content by means of equation (19) suffers 
especially near the point where This is easily seen by considering the 

equation 


* a„ Br_ + 


0 n c cos ()> 


o ) ^ ^ ^ j ^c *^0 *^0 ^ ^ '' 


If the orbit is circular, then the expression for (l/z^,) J N dz in terms of is 


( 20 ) 



b* 


1 c ^ 

sin <>Q cos <>Q J ax 


z. dz . 


( 21 ) 
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Ideally, in the abaence of KmdientM, M -» 0 as 0, niid it 1 h poHaible to calculate 
dS near the vertical by numerical application of the I'lioopital rule. In retility, 
however, vi/e obtain the value 



“ o * c **"*0 


( 22 ) 


which may include a air.nble horizontal -((radient eifect. VMten the Rutellitc ia directly 
over the station (<>q” 0), the differential duppler Hhlft will l>e approxlmatelv 
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or simply a function of the horizontal gi.; iient. Again, this is on idealization, since 

the higher order effects neglected in deriving the abo\e equations will at times add 

• 

a few percent to the theoretical value of <>«&. One example is that of the overhead 
(<))q ■ 0) case. The doppler shift of each frequency depends not on the straight-line 
relative velocity between the antenna centers but on the component in the direction of 
the wave normal at the emitter (Al'pert, 19U:i). The difference becomes important 
near the vertical and will l3e manifested ns a small displacement from the vertical of 
the direction where 6* » 0. We mention these effects primarily liecause of their 
Importance in the idealized simulation of f>4> and in the subsequent reccnerv of 
ionospheric parameters. The actual oiiser.’ntions will have components due to out-of- 
plane motion, orbital eccentricity, and nonstationarity of the ionosphere (BN/dt ^ 0) 
predominating over the small second-order omissions in ^he theory. 


In order to demonstrate what we c:m expect in the actual data reduction, we have 
simulated the DM-to-ground differential doppler records for three widely varying 
ionospheric conditions. Equation (8) is simplified by a circular orbit of the satellite 
(z^ = 0) with the ground station always in the orbital plane (y^ = 0). Furthermore, 
the effect of nonstationarity |9N/Bt) is excluded. For the simulatiffli, we chose to 
locate the ground station at positions corresponding to 0430, 0930, and 1230 LT. At 
0430 LT, the electron density is lowest, with correspondingly small gradients, and is 
the least interesting case for our purposes. The small doppler shifts simulated are 


22 



subject to a larger relative measurement error, and computer time limitations on the 
ray-tracing accuracy prevented ue from making any clear distinction of the small 
gradient effects. Wo will therefore ctmcentrate on the 0930 and 1230 LT cases, wlicro 
the horizontal gradients are clearly significant. 


One way to examine the effect of neglecting 8N/0x in toe Inversion Is to compare 
the total simulated with the theoretical value, given by 



ir 

I c cos 



(ii4) 


where the integral is computed directly from the ionospheric model. We shall call 

6^ the "columnar" dopplor shift, to denote Its dependence only on the average columnar 
c 

electron content between the satellite and the ground station. Figures 8a, 8b, and 8c 
show the relationship between b4 and 6$^, for various elevation angles with the station 
at 04?0, 0930, and 1230 LT, The elevation E is measured positively from toe horizon 
over which the satellite approaches toe station. 


The inversicm of our expected 6h to columnar electron content without 

accounting in any way for the l9S/3x) term is perilous indeed. In Figure 9a, we have 

plotted the ejected actual columnar electron content for the 1230 LT DM-to-ground case 

(cur\'e A), and it should be compared with curve B, which gives toe columnar electron con- 

f 

tent obtained by simply inverting equation (19) for N dS. There is a sizable error 
throughout the pass and a noticeable divergence near the vertical, owing to the fact that 
6f> aimroaches not zero, but the value given by equation (23), By applying toe 1' Hospital 
rule at 4>q = 0, we get a value of (Sj^j + laN/axJ^; the corresponding apparent columnar 
electron content is identified by Bj^ in Figure 9a. 

At 0930 LT, the horizontal gradient is much larger, as evidenced by the value at 
221 km (Figure 3), and the effect on curve B (Figure 9b) is severe. 

An Immediate correction can be introduced assuming J(8N/0x) z dz to be con- 
stant throughout the satellite pass if 64^ = ((jj^ = 0) can be detected. Within the first- 

order planar- approximation theory. 
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Figure 8a. Expected differential dealer shift in the DM -to-ground path at 0430 LT. 
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Figrure 8b. Expected differential doppler shift in the DM -to-ground path at 0930 LT. 
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Fi^re 8c. Expected differential doppler shift in the DM -to-ground path at 1230 LT. 


26 




0930 LT 
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Figure 9b. Columnar electron content in the DM -to-ground patli at 0930 LT 


( 25 ) 


• * H 


If the integral in equation (25) is a function of only, then its cuntriiiutian to 64> is sim- 
ply fxf^/cos which can ije subtracted Ijefore so!*.ing for the columnar electron content. 
By interpolation about 0, we found a value of 6*^^ ■ +38 mHz for the 0»30 LT pass 
and a smaller negative value of -8.7 mHz for 1230 L7. Applying equation (25) and 
using the integral calculated from our model ionosphere (Grossi an** Gay, i975, Api>en- 
dix A), we found the corresponding theoratical values of 64^ to l)e -7 jm.i *39.3 mHz. 
Hence, if the 6«f> record is sufficiently smoothed near 4)^=0, we can obtain (4>(0) rather 
accurately in the presence of a lax'ge horizontal gradient. The corrected columnar electron 
content is represented by curve C in Figures 9a and 9b, This correction succeeds in 
removing a large part of the error near 4)^= 0, but it is less effecti\e at larger .mgles 
since it does not take into account the longitudinal variation of z dz, an 

effect that can be seen at 0930 LT, where the gradient contrilnjtion on the left (E < 90*) 
is much smaller than that on the right. Hence, cun-e C reflects an overcorrection to 
the columnar electron content on tlie left-h'Uid side. In Figure 9a, the coriection is 
actually in the wrong direction for E > 90' but reduces the error considerably on the 
left. 


4,2 Inversion by Use of Simply Modeled Horizontal Gradients 

A suitable ionospheric electron-density mo^lel incorporating knowledge of horizontal 
gradients can i)e emplo>'ed to estimate the effects of dN/dx ;md 9N/9y on the obser\e<i 
64>, The layer models of Chapm:m (1931) have a wide application in this l>pe of prfx;e- 
durc, and several variations can tie considered. The selection of a particular model 
should be affected by the most current reliable obsenatlons and predictions of 
ionospheric electron density for the time of the ASTP mission. 

In our simulation of data reduction, we have chosen the simplest of models, a 
Hr jar variation of 9N/9x with altitude. Since the numijer of free electrons (and 9N fix) 
diminishes to nearly zero at heights of 70 to 80 km, we chose a linear profile, as 
depicted in Figure 10, where 9N/9x is zero i>elow 70 km and takes the value of (9N/9x)| 
given by the model used in the simulation (Grossi and Gay, 1975, AppendLx A). 
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Figure 10. Linear gratlient profile. 
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The results obtained by incorpomting these linear-nuxiel Krndients into the cal- 
culation of the total electron content ars given in Figures llo and lib. The curves • 
labeled D represent the columnar electron content obtained from fa* after corrections 
have been made for the effect of iB^i/dxl computed from the aix)ve model. In the 1230 LT 
case (Figure 11a), curve D would fall Ijetween curves H and C in both quadrants in 
Figure 9a. The correctivjn based on the a priori linear model is then somewhat less 
effective than the correction applied by using fa*^. 

At 0930 LT (Figure 11b), curve I) if plotted in Figure 9b would fall below curve C 
during the satellite's approach to the ground station, and in the receding quadrant, it 
would fall between C and the actual columnar-electron-content curve (A). It is ap{)nrent 
that the change in (oN/Oxl with elevation (luring a pass is sigrtlficant imd that a mf)re 
sophisticated model is needed. 

4.3 Application of Spacecraft -to-9pacecraft Horizontal Gradient Measurements 

Instead of relying completely on molels of ON /fax, we are able, in the ASTP experi- 
ment, to utilize the actual measurements of the horizontal gradient in the DM/CSM path 
in the data reduction. It is well knowTi that the ionospliere is not ideally smooth over the 
characteristic observing period of 2 to 3 min required for a satellite pass 
(Al'pert and Sinel'nikov, 19G5). Irregular short-|)eriod variations in fa«t somewhat 
obscure the smot'dh average electron-density characteristics in the region of the 
satelllte-to-ground path. From DM/CSM measurements of fa«t, we can olitaln know'ledge 
of the irregularity of N along the orbit, limited only by the 10-sec (77-km) resolution 
provided by the on-board dv^pler records and the 300-km distaivce l)etwtHjn the space- 
craft. Since the ASl’P orbit is nearly circular, the Irregular variations due to 
will be small. The possibility also exists that the DM/CSM ionosph<*ric doppler varia- 
tions can lie used to identifv clearly the effect of N z :md enable us to remove that 
term from consideration in the DM-to-ground doppler shift. 

In computing ionospheric parameters via equation (8) from observations of fa* in 
the spacecraft -to-ground case, knowledge of the magnitude of the horizontal gradient 
is essential to an accurate description of the ionospheric electron content l)Ctvveen the 
transmitter and the receiver. Although the single measurement of (iS at the 
orbital height gives no direct information about the vertical profile of 9N/ax, it may 
yet be more important since the integral (9?5/ax| is vzeighted by the height z. 
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Figure 11a. Continuation of Figure 9a. 



0930 L i‘ 
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Figure lib. Continuation of Figure 9b. 




To (lexnonFtrnte what we expect to achieve from the aloft gradient measurement 
In terms of Inversion accuracy, we have again usei! a linear profile of HN/ax iietween 
221- and 70-kin altitude, with the upper values lielng determine<l by the oi)served 
3N/8X Ijetween the DM awl the CSM. That expression Is 

yN I » '/. 2 70 km , 

aT * ^ 0 , z s 70 km , 

where x Is the <listance lietxv’eon the subsatellite point an<l the ground station, anti 

2 

Zq “ 70 km. The slope a(x) » 9 N/9z9x Is a quadratic determlne<l from the values 
taken above the ground station at points along the orbit corresponding to 15* of longi- 
tude either side of the station. 

Referring once again to Figures 11a and 11b, we see the results In terms of columnar 
electron content displayed by the curves lalicled E. At 1230 LT, when the satellite Is 
approaching the station, the E cune is almost coincident with C as given in Figure 9a 
and Is an Improvement o\er D. In the receding quadr:mt, however, the E and I) cur es 
are displaced by a rather constant amount from the actual columnar electron content, 
suggesting that the vertlc;il profile of 9N/9x is more poorly modeled than the horizontal 
variability. 

For 0930 LT (Figure 11b), the E cune is nearly the same as the D cur e on the 
left (satellite approaching) and diverges slightly near the vertical. On the right, the 
E values of the columnar electron content are l)etter than those of the 1) cune, and again 
there is still some divergence near the vertical. It is interesting to note that in the 
left-hawl quadrant, the cunes C, D, and E (Figures 9b and 11b) are successively 
rerno\ed from the actual columnar electron content. A, while on tlie right, the opposite 
is true. That is, the sequence of corrections (C, D, E) we have applied causes the 
solution to diverge from the real value in one quadrant and to converge in the other 
quadrant. 
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5. CONCLUSIONS 


Tto analysis of the expected DM-to-CSM differential dq^lor shifts Indicates that 
a measurement accuracy of 1 to l(yj in the value of the horizontal electron-density 
gradient at 221-km altitude can he achieved, 'llio accuracy increases with the gradient 
amplitude and is sufficient in each case cited to account adequately for the ionospheric 
frequency shift in the DM/CSM link, which must be removed for the gravitational experi- 
ment, Possible exceptions {where the accuracy of 9N/ax is insufficient) are at the 
nighttime minimum of electron density and in the presence of large temporal variations 
in N along the radio path. 

In the case of ground-based obsen’atlons near the orbital plane, we have Identified 
some possible effects of the horizontal gradient 0N/bx on the received differential 
dc^jpler shift and have successfully applied a simple method of compensating for the 
presence of horizontal gradients in the solution for space -to- ground columnar electrcai 
content, m the two cases considered, it was possible to reduce the gradient-associated 
error in the columnar electron content by approximately one-half by introducing a 
simple, linear, vertical-profile model of SN/ax. A more aj^ropriate gradient model 
with more horizontal and vertical details should result in an overall accuracy of 5 to 
10% in the columnar electron content. In our simulaticm, the horizontal -gradient 
measurement at the ASTP altitude employed in conjunction with a linear model of the 
gradient profile gives essentially the same correction as does the a priori linear 
model. The rms difference between the simulated columnar electron content and that 
given by cur\'e E in figures 9 and 11 is G, 5-y in tlie 1230 LT case and 7, 5% for the 
0930 LT case. 

The results of these numerictil experiments, performed imder assumed idealized 
conditions, cannot yet be applied broadly without further consideration of the effects 
we have neglected, especially at satellite altitudes different from the ASTP orbital 
height. We hope to perform more simulations tliat will include the earth's magnetic 
field, the time variability of N, and the effects of out-of-plane motiem (y^ t 0) if time 
permits. 
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